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Abstract — In this work, an analytical model to calculate the 
channel potential and current-voltage characteristics in a Sym- 
metric tunneling Field-Effect-Transistor (SymFET) is presented. 
The current in a SymFET flows by tunneling from an n-type 
graphene layer to a p-type graphene layer. A large current peak 
occurs when the Dirac points are aligned at a particular drain- 
to-source bias Vds- Our model shows that the current of the 
SymFET is very weakly dependent on temperature. The resonant 
current peak is controlled by chemical doping and applied gate 
bias. The on/off ratio increases with graphene coherence length 
and doping. The symmetric resonant peak is a good candidate 
for high-speed analog applications, and can enable digital logic 
similar to the BiSFET. Our analytical model also offers the 
benefit of permitting simple analysis of features such as the 
full-width-at-half-maximum (FWHM) of the resonant peak and 
higher order harmonics of the nonlinear current. The SymFET 
takes advantage of the perfect symmetry of the bandstructure 
of 2D graphene, a feature that is not present in conventional 
semiconductors. 



I. Introduction 

Graphene is an atomically thin two dimensional (2D) crystal 
0. Due to the high mobility of carriers in it, their linear 
dispersion and perfect 2D confinement, graphene is being 
considered as a channel material for future electronic devices. 
However, some challenges such as opening a finite bandgap 
for digital applications still remain. 

A distinguishing feature of graphene is its symmetric elec- 
tronic bandstructure. The valence band is a perfect mirror 
image of the conduction band about the Dirac point. Such 
symmetry carries over to gapped 2D crystals such as hexag- 
onal Boron Nitiide (h-BN), and less so to transition metal 
dichalcogenides (such as Molybdenum Disulfide, MoS2)|2|. 
This unique bandstructure symmetry has not been sufficiently 
exploited for active device applications till date. 

Most studies of graphene based devices have focused on car- 
rier transport in the 2D plane of the crystal. Recently, however, 
carrier transport out of the plane, i.e., vertical to the graphene 
sheet, has received increased attention. These studies of out-of- 
plane charge transport in 2D crystals have been motivated by 
the proposal of the bilayer pseudo-spin FET (BiSFET) in 2009 
0. The BiSFET exploits the fact that two graphene layers can 
be placed in close proximity and if populated by elections and 
holes, the strong Coulomb attraction between them can lead 
to exciton formation. Excitons are bosonic quasiparticles, and 
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can undergo condensation below a certain critical temperature. 
Since the Fermi degeneracy in graphene is tunable over a 
large energy window, the critical temperature for the excitonic 
condensate has been calculated to be higher than room temper- 
ature. The formation of the excitonic condensate is expected 
to lead to a macroscopic tunneling current between the layers. 
Similar behavior has been observed at low temperatures and 
at high magnetic fields in coupled AlGaAs/GaAs electron- 
hole bilayers 01,0. The BiSFET thus has the potential to 
realize many-body excitonic tunneling phenomena at room 
temperature. The power dissipation in computating using the 
functionality of BiSFET is predicted to be many orders lower 
than conventional CMOS switching. 

Stacking of different 2D crystals leads to a novel class 
of heterostructures [6 |. For example, stacking graphene with 
BN results in a smooth surface, since BN shares the same 
hexagonal lattice structure with graphene. The absence of out- 
of-plane covalent bonds implies that strain effects in similar 
lattice mismatched heterostructures based on sp 3 -bonded 3D 
crystals are much reduced, or perhaps even eliminated. At 
low carrier concentration n ~ 10 11 cm -2 , the device exhibits 
mobilities of the order of 100,000 cm 2 /V.s at room temper- 
ature, which is much higher than for graphene on Si02 or 
SiC [7|. Electron tiansport out of the plane of graphene has 
also started receiving attention in experiments. For example, a 
graphene/BN/graphene sandwich heterostructure was recently 
reported [8|, and interlayer election transport was measured. 
This stacked graphene/BN/graphene heterostructure showed 
a room temperature switching ratio of 50 and 10,000 for a 
similar graphene/MoS?/graphene structure [8 |. Another recent 
report [9| showed that electrons can be moved from 2D 
graphene into and out of Silicon to form a variable barrier- 
height device (the Barristor). 

As charge transport out of the plane of graphene receives 
increasing attention, a pertinent question emerges upon careful 
analysis of the proposed BiSFET device. What is the expected 
behavior of a similar device stiucture consisting of a graphene- 
insulator-graphene (GIG) p-n junction heterostructure, but in 
the absence of the many-body excitonic condensate? Similar 
devices already exist in III-V resonant tunneling diodes, where 
single-particle tunneling itself leads to a number of interesting 
and useful quantum phenomena that persist at room tempera- 
ture. Negative differential resistance is one such effect. Single- 
particle tunneling current transport has found enhanced atten- 
tion recently in homoj unction and he teroj unction Tunneling 
Field-Effect Transistors (TFETs). It has been measured across 
various semiconductor heterostructures at room temperature, 
highlighting its robustness iflOl . 

One major novel feature of graphene is the perfect symmetry 
of the bandstructure, which can lead to enhanced functionality. 
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Fig. 1 : (a) Sketch of the SymFET and the energy-band diagrams for a doped graphene/insulator/graphene junction, at voltages 
of (b) qVos < 2AE, (c) qVos > 2AE, and (d) qVos - 2AE. A qualitative current- voltage I—Vds characteristic is shown in (e). 
The inset in (a) shows the symbol defined for the SymFET. 



Motivated by the above question, we recently calculated the 
single-particle interlayer tunneling current-voltage curves ex- 
plicitly for finite area two-terminal GIG heterostructures ifTTI . 
The general finding was that at most interlayer bias voltages, 
energy and momentum conservations force a small tunneling 
current to flow at one particular energy halfway between the 
two Dirac points. However, at a particular interlayer voltage 
when the Dirac points of the p- and n-type graphene layers 
align, a very large interlayer tunneling current flows. This is 
because energy and momentum are conserved in this process 
for all electron energies between the quasi-Fermi levels of the 
n- and p-type graphene layers. The I—V curve is dominated 
by a Dirac-delta function-like peak at the critical interlayer 
voltage, and smaller currents at all other voltages. Our explicit 
calculation of the tunneling current also showed that the effect 
is highly robust to temperature, but less robust to rotational 
misalignment of the two graphene layers. This surprising, yet 
conceptually simple behavior of the 2-terminal GIG device 
leads naturally to the question: how will a transistor geometry 
with the single-particle GIG tunneling junction as its channel 
behave? 

In this work, we extend the detailed physical model of the 2- 
terminal GIG device described in IfTTI to a Symmetric graphene 
tunneling Field Effect Transistor, which we call the "SymFET" 
since its unique characteristics derive from the symmetry of 
the bandstructure. We derive analytical expressions for the 
channel potential and current-voltage characteristics of the 
SymFET. Possible logic and high frequency applications are 
also discussed. 

II. Device Model 
We assume a symmetric device structure as shown in 
Fig. Q2 a )- Two graphene layers are separated by an insulator, 
and this GIG structure is sandwiched between a top and bottom 
gate. Ohmic contacts are formed to the two graphene layers 
individually representing the source (S) and the drain (D). The 



top and bottom gate voltages Vjg, Vbg control the quasi-Fermi 
levels fi„ and fi p in the top and bottom layers of graphene. The 
gate insulator thicknesses of both gates are assumed to be the 
same. The quasi-Fermi level is AE above the Dirac point in 
the n-type graphene layer and below the Dirac point in the 
p-type graphene layer. This is indicated in Fig. Q2b), (c), and 
(d). The top and back gates are symmetric Vjg = ~Vbg< an d 
the drain-source voltage is Vds — Vd - Vs . The insert in Fig. Q] 
(a) shows a proposed device symbol for the SymFET. 

As shown in Fig. |TJb) and (c), under S/D biases when the 
two Dirac points are misaligned, only a single energy (and 
lateral k-momentum ring) in the Dirac cone meets the require- 
ment of simultaneous energy and momentum conservation, 
and thus the tunneling current is small. At Vds = 2AE/q, 
however, the two Dirac points align, and electrons at all 
energies between the quasi-Fermi levels satisfy energy and 
momentum conservation. A large tunneling current is thus ex- 
pected; a resonant current peak originating from the perfectly 
symmetric bandstructure of the graphene layers should result. 
Since graphene is not a metal, a part of the applied voltage 
will drop in the graphene layer itself. This effect of the finite 
density of states is captured in the quantum capacitance of 
graphene and is included in the model. This is critical since 
the gate capacitance C g = e g /t g and the tunneling capacitance 
C t = £t/tt are large and can reach the quantum capacitance 
limit easily lfl2l . For simplicity of the calculation we use the 
T — > OK limit for quantum capacitance: 



c a = 2 - 



\AE\ 



(1) 



n(hv F /q) 2 ' 

where q is the single electron charge, Vf is the Fermi velocity 
in graphene, ft is the reduced Planck constant. It can be verified 
that this approximation is a good one at room temperature. 

The source and drain electrodes are assumed to be perfect 
ohmic contacts for simplicity. In practice, the contact resis- 
tance will also force a voltage drop, and can be added on top of 
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the intrinsic model. The interlayer tunneling current calculated 
in IfTTI depends on the interlayer voltage difference. Therefore, 
to find the behavior of the 4-terminal SymFET, potentials of 
the two graphene layers (p c hi and <p c ta need to be identified 
as a function of the gate and drain/source biases (the channel 
potentials of graphene and Fermi levels are referenced to the 
aligned channel potentials at flat band.). To do so, we invoke 
the charge neutrality condition: 



electron mass, d is the normalization constant for z-component 
wavefunction in graphene, uyi is a constant of order unity, and 
A = L 2 , with L being the coherence length of graphene (size 
of ordered areas in graphene film). In this work, we assume 
the graphene size is A with coherence in all area. 

The resonant current is a perfect Dirac-delta function at 
qVos - 2AE for infinitely wide graphene sheets. For finite 
widths of length L, it is broadened to IfTTI : 



( + Vtg)C s + ( )C, + ( )C q /2 + qN = 0, 

q q q q q 

(2) 

( + Vbg)Q + ( )C t + ( )C g /2-qN = 0, 

g q q q q 

(3) 



-Vbg = Vg, AT = AE 2 opjng /n(tiv F ) 2 is the chem- 



where VYg 

ical doping concentration, and we assume the work functions 
of the metals are matched with the undoped graphene sheets 
which gives the flat-band conditions at zero gate bias. The 
factors 1/2 in the third terms of equations (fj) and (01 are 
due to the linear dependence on AE in the graphene quantum 
capacitance, which is a differential capacitance. 

Then we take equation (f2]) minus equation (|3J, use relation- 
ships: qV DS = 2AE + tp ch \ - <p ch2 , qV DS = Hn -fJ. p , Hn - <p c hl = 
AE, and ip cn % - Hp — AE, we can form a quadratic equation 
with the only unknown parameter AE: 



AE AE 2qAE 2 

(V DS -2— + 2V G )C H + 2(V DS - 2— )Q - - - + 2qN = 0, 
q q nyhvpY 

(4) 



and the solution is 



AE(V G ,V DS ): 



(2C t + C g MHv F /q) 2 



(2C t + C g ) 2 7T 2 (Hv F /q) 4 



n(hv F y 
2q 
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[(Yds + 2V G )C g + 2C,V DS + 2qN] 



(5) 

The electrostatic model used here is based on a ID approxi- 
mation, ignoring the intra-graphene layer potential distribution 
and current flows. A more rigorous treatment requires the 
solution of the 2D Poisson equation, which is suggested for 
future work. 

The analytical expression for the interlayer tunneling current 
at zero temperature was derived in [1 lj. When the Dirac points 
in the two graphene layers are misaligned, the nonresonant 
tunneling current is: 



2A£ 

= Gi( Vds), (0 < qV DS < 2AE), 



2AE 



(6) 



/ = Gi(Vds ), (qV DS > 2AE or qV DS < 0), (7) 



q 



<f-AMAl e *"\2 



where the prefactor conductance G\ = ^-( —jj^fc- — Y, k is a 
decay constant for the tunneling current in barrier, m is the free 



1.6 LAE 2 (2u A n +u\ 2 ) 



V2lr u\ 2 qhv F 



A qVDs-2AE 2 

exp[-— ( ) ], (8) 

47T nv F 



where u\\ is a constant of order unity, similar to u\2. The 
total current is the combination of Equation (|6]l or Equation 
(|7]i summed up together with Equation ©. 

Equations ([6]) - (O were derived assuming zero temperature. 
At finite temperature, to capture the thermal occupation of 
states, the current needs to be calculated by including the 
Fermi-Dirac distributions in the integral over all states IfTTII . 
Equations (O and ( fTOb are the finite temperature expressions 
corresponding to Equations (|6) and ©, derived by a direct 
extension of the theory in our previous work IfTTI . 



I = G 



4H 2 v 2 F 



J {k[f{E n , k -Mn, T) - f(E pJl -fX p ,T)] 



: 6(2AE - qV DS - 2hv F k)\dk 



(9a) 



2A£ 



: Gl( V DS )[f(-qV DS /2,T)-f{qV DS /2,T)] (9b) 

q 

2AE qV DS 

■ Gi( V DS )tanh(^), (qV DS < 2AE) (9c) 

q 4k B T 



4h z v 



2, ,2 



J [k[f(E nM -fi„,T)-f(E p , k -h p , T)] 



'- 6{qV D s ~ 2AE - 2hv F k))dk 



(10a) 



■Gy(V DS - 



2AE 



■)[f(-qV DS /2,T)-f(qV DS /2,T)] (10b) 



q 

2AE qV DS 

-)tanh(^-^), (qV DS >2AE) 



(10c) 



q 4k B T' 

where subscripts n and p refer to the n-type (top) and p- 
type (bottom) electrodes, respectively, and with the Fermi- 
Dirac distribution for arguments SE and T given by f{5E, T) = 
1/[1 + exp(6E /ksT)]. Following Ref. IfTTI we have for Eq. 
(|9af >, with < qV DS < 2AE, that E n<k -fj, n = fiv F k — AE and 
Epji-fip = -hv F k + AE, hence yielding Eqs. d9bl ) and d9c| ). 
Those equations also hold for qVos < 0, with the difference 
between the Fermi-Dirac occupation factors changing sign. For 
Eq. ( I lOab . with qVos > 2AE, we have E„^ -pin - —%v F k- AE 
and Epfc—Hp = Hv F k + AE, yielding Eqs. ( II Obi ) and 111 Oct , 

A finite-temperature correction to Eq. (|8) is more 
complicated, since that equation itself is a significant 
approximation IfTTI . Nevertheless, the dominant term in its 
temperature dependence can be recognized as the increased 
number of states available for tunneling, in the resonant 
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Fig. 2: Id ws.Vds characteristics showing scaling with (a) 
the tunneling insulator thickness and (b) the gate insulator 
thickness (an arbitrary k value is chosen here to give a clear 
illustration of current density scaling). 



0.5nm 



O.OeV : 




(a) 



(b) 



Fig. 3: The contour plot of the complete bias space for (a) 
chemical doping AEdoping = O.leV and (b) no chemical doping. 



situation of Fig. 1(d) with qVos - 2AE. For this resonant 
situation the Dirac points of the two electrodes are aligned, 
at energy E^. This number of states is given simply by 



-t-oe 

■/■ 



N S (T)= I p(E-E d )[f(E-E d -AE,T)-f(E-E d +AE,T)]dE, 

(11) 

where p(E) = 2\E\/[7r(Hvf) 2 ] is the density-of-states (DOS) per 
unit area of graphene. Our correction to Eq. (|8) is then made by 
multiplying it by N s (T)/N s (0), with N s (0) = 2AE 2 /[n(tiv F ) 2 ]. 
Expressing the result in terms of Fermi-Dirac integrals of order 
1, 3i, we have 



N S (T) 2(k B T) 2 



AE 



AE 



N s (0) 
where 3 1 (x) = 



AE 2 



[3i(— )-3i(-— )] 

KbT IcbT 



(12) 



f/[l +exp(t- x)]dt. 



Equations (0 and ( fTOt , as well as the product of Eqs. ([8} 
and (fT2l . are valid only in the limit of large L, although the 
L value at which the approximations break down is different 
for the two cases. In the former case, the dominant term in 
the slope of I(Vds) at Vds - is qVos l^kaT, arising from 
the tanh terms, and this slope becomes very large for small T. 
However, the exact solution for the current [11] yields a slope 
that is limited by the L value. We find that simply multiplying 
Eqs. (0 and (fT0T > by a factor of tanh(LqVos /ttHvf) yields a 
slope at that agrees very well with the exact solution, and it 
does not significantly affect the I(Vds) curve elsewhere. For 
the product of Eqs. (O and (flZb . the current from Eq. (|8) is 
nonzero at Vds=Q > with this discrepancy being significant 
only for sufficiently small L values. In this case we find that 
multiplying the product of Eqs. (© and (fT2l by a factor of 
tanh(LqVos /2nhvf) solves this problem of the nonzero current 
at Vds = , and it also produces a slope at Vds = that agrees 
fairly well with the exact solution. Hence, our final formula 
for the total current is given by 



2AE 2AE qV D s LqV D s 

I = G^Vds ~ —)sgn(V D s ~ — )tanh(^)tanh(4-^) 



+ 



q q 
f 1.6 ^ LAE 2 (2u 4 n + u\ 2 ) 



nhvf 



I V2^ 

NAT) 
NA0) 



u 4 i2 qhvp 
tanh(— — ) 



A qVos-2AE 2 

exp[-— ( — ri 

4k nvf 



(13) 

where sgniV^s -2AE/q) equals 1 for Vds > 2AE/q, for 
V DS = 2AE/q , or -1 for V DS < 2AE/q. 

The single particle tunneling model used in this work cap- 
tures the relevant physics e.g. wave function overlap (detailed 
derivation in Ref. ifTTI '). even in the case of strong inter-layer 
interactions at small insulator thicknesses. The only approxi- 
mation we made is that the electronic structure (band structure) 
of the GIG system is the same as two non-interacting graphene 
sheets. We have justified the approximation as follows: any 
modification to the band structure due to the interaction will 
be near the Dirac energy. At sufficiently high doping AE, the 
error is negligible. 

III. Results and Discussions 

We present typical I -V characteristics of the SymFET at 
room temperature first. The comparison of T — 300 K and 
T — K is discussed later. The values of the decay constant 
k can be calculated from the complex band structure inside 
the bandgap of the insulator based within the effective mass 
approximation HQ. Here, we use an estimated value k — 17 
nm -1 , following footnote 14 of IfTTI . The chemical doping 
level is set to be AEdoping = 0.1 eV. A finite coherence length 
L — 100 nm is assumed; the effect of this parameter on the 
on/off ratio for the device is discussed below. When the Sym- 
FET scales down to 50nm or less, momentum conservation 
does not scale well with device size. Scaling limits of SymFET 
will be a subject of future study; the quantum confinement and 
quantized transverse momentum in graphene nanoribbon also 
need to be considered near scaling limits, as discussed in Ref. 

ma. 

The resultant Id -Vds characteristics with varying insulator 
thickness are shown in Fig. |2] The tunneling insulator thick- 
ness t, is similar as the tunneling barrier thickness in double 
quantum well heterostrucures [5 |. As t t increases the resonant 
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Fig. 4: (a) vs. V/^ curves with different Vg, and (b) vs. 
Vq at different V/^. 



peak current decreases, as expected. The gate insulator can 
be a high-k material similar to that employed in Si CMOS 
technology. In addition, 2D materials such as BN might be 
a better choice to reduce the interface trap density since the 
dangling bond can be reduced. The measured breakdown field 
is as high as 7.94 MV/cm for BN 03). Thinner t g offers 
better gate control and higher gate induced doping. When t g 
decreases, AE becomes larger at same gate bias. The resonant 
peak moves to a higher bias and the peak current increases. For 
the simulation results shown next, we fix the gate capacitance 
with t g = 20 nm and dielectric contact e g = 9, and the tunneling 
insulator thickness t t = 0.5 nm, and dielectric contact e, = 9. 
(BN might have an even lower dielectric constant €bn = 3.5 

ESI). 

In Fig. [3] the entire bias phase space of the I—V charac- 
teristics are shown. When the graphene sheets are chemically 
doped (i.e., AEdoping =0.1 eV), Vds drives the carriers tunnel- 
ing between p- and n-type of graphene and the resonant peak 
exists for Vg = V . With non-zero Vg, gate electrostatic 
doping further increase AE and resonant peak shift to higher 
Vds ■ But the current is quite small in the non-resonant region. 

If we define the peak current as the /„„ and the current close 
to Vds ~ as I ff, from Equation © and ([H) the effective 
on/off ratio is: 



Ion 
I off 



0.8 LAE 



(14) 



V2tt hv F 

The effective on/off ratio is independent of temperature 
(ignoring the slight difference due to the Fermi tail), and 
increases with the doping AE and the graphene size L. For 
L = 100 nm the ratio is -100 and -1000 for L- I /urn. We 
also point out that, the on/off ratio is not necessarily a figure of 
merit for the device (since the device might be employed for 
analog applications where high modulation is not required). In 
a similar vein, the resonant peak is symmetric in voltage, and 
represents a rather strong negative differential resistance. The 
peak-to-valley-current-ratio (PVCR) of this NDR is identical 
to the on/off ratio defined above. 

In Fig. [3](b), we assume that AE doping = eV. At V G = V, 
AE is non-zero since drain bias also induces doping, similar 
as the case of GIG junction in 11 11 . However, AE is small 
and the resonant peak is small at low Vg- When Vg increases, 
electrostatic doping induces an appreciable resonant current, 



which further increases at higher Vg- 

The Id -Vds characteristics at fixed Vg is shown in Fig. [4] 

(a) . The resonant behavior shows clear on and off states 
without a saturation region. Because we assume a chemical 
doping of graphene, the SymFET with a resonant current peak 
can operate at Vg = 0V. As mentioned above, the gate will 
induce electrostatic doping in the graphene layer. With larger 
Vg, AE increases, the resonant condition qVos = 2A£ occurs 
at larger drain bias, and the resonant current peak moves to the 
right. Higher Vg induces more doping and thus large on-state 
current. In the 2-terminal GIG device, the resonant current 
peak is proportional to the coherence length L, and the width 
is proportional to 1/L IfTTI . In the gated SymFET, since the 
gate bias electrostatically dopes the graphene, it offers the 
additional flexibility to adjust the on and off states. In Fig. [4] 

(b) , Id - Vg curves are shown with a strong non-linear and 
resonant behavior, but with wider peaks. When Vg is small 
and outside the resonant peak, the transconductance is small, 
but is large in the peak condition. 

Different from the sub-threshold region of MOSFETs, we 
refer to the 'off' state away from the resonant peak in the 
SymFET as the non-resonant region. Although the graphene 
sheets are doped, according to conservation laws, the current 
near V G = V is low (V D s = 0.7 V, solid line in Fig.E(b)). The 
gate voltage modulates the doping potential AE in SymFET. 
By Eq. Q, at fixed Vds, doping potential AE is a sub-linear 
function of Vg- The current Ids(Vg) at non-resonant region 
roughly follows the same hyperbolic form as the first term of 
Ids (Vds) in Eq. (13]): 



/octanh( ^)tanh(^). 



(15) 



Hyperbolic functions may not offer a sharp subthreshold- 
swing (SS), but the SymFET is more attractive for analog 
applications, where the steep SS is not necessary. 

The resonant current peak follows a normal distribu- 
tion function / ~ exp(- teYMz£Ml )^ with Full-width-at-half- 



maximum, FWHM = 2.3548cr, where <x = 



for the 



resonant peak in SymFET. In fact, AE is not a constant but 
also dependent on Vds and Vg, thus the FWHM shown here 
is only an approximation. The reason why a narrower FWHM 
occurs for Vds compared with Vg in Fig. [4] relates with the 
large gate insulator thickness (weaker gate control). Smaller 
FWHM of Vg can be achieved with thinner gate insulator (not 
shown here). 

Because tunneling is the main current transport mechanism, 
the Id - Vds curve is quite insensitive to temperature as shown 
in Fig. [5] But since the Fermi Dirac distribution smears out 
state occupancy at a finite temperature, slight differences can 
still be observed between T = 300K and T = OK. At low V D s, 
the transport energy window (between the quasi-Fermi levels 
ju„ and /dp) is small (Fig. [5] (b)). Then, the Fermi distribution 
smearing reduces the carrier density at higher temperature 
and the current decreases. The increase of the resonant peak 
current at room temperature is because the Fermi distribution 
tail extends to higher energy with more states. When the 
Dirac points are aligned, states at all energies conserve lateral 
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Fig. 5: (a) The comparison of Id~Vds characteristics at T — 
300K and T = OK. (b) the sketch of the band diagram and 
Fermi Dirac distribution at T = 300K and T = OK. The inset 
of (a) shows the current density near Vds - 
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Fig. 6: The nonlinear resonant current is highly symmetric. 
When a DC voltage Vdsp biased at the current peak is 
superimposed with an ac signal Vds(t), the frequency of the 
output current i(t) will be doubled. 



momentum upon tunneling, and thus are allowed. 

We note that the device current has symmetric resonances in 
both the Id~ Vds an d Id — Vg scans. This is quite unlike what 
happens in a single-layer graphene FET where the 'ambipolar' 
nature manifests itself primarily in the gate bias sweep lfl6l 
irTTI . The nonlinear symmetric resonant Id - Vds behavior can 
be used for purpose of frequency multiplication (Fig. |6). If a 
DC voltage bias at the current peak Vosp is superimposed 
with an ac signal, the frequency of the output current will 
be doubled. We can use equation dS) to calculate the output 
ac current signal. Assuming that Vds — Vdsp + Vds elcot > the 
oscillatory part of the current is: 



1.6 LAE 2 (2u 4 u +u 4 2 ) 



V2^ 



u. 2 qhvf 



exp(- — [ 



A r v ds e" 



-?)■ (16) 



To find out the higher order harmonics, we ignore the constant 
prefactor: 



/aexp(- — [— Y). 

4tt nvF 

This expression can be further extended as: 



(17) 



C 2 C 3 
I a 1 - Ci exp(2;wf) + -j- exp(4/wf) - -i exp(8;wf) + • ■ • . 



(18) 

jV 2 , . In equation dT8l i. only even higher 



where C\ = - 
harmonics occur. 

The SymFET is expected to be intrinsically fast since it 
relies entirely on tunneling. The extrinsic performance with 
parasitics can be analyzed same as for any high speed device 
and will not be covered in this paper. High-frequency digital 
operation and a host of analog applications such as frequency 
multiplication are thus possible by exploiting the symmetry of 
the bandstructure of 2D graphene. 

As explained in ifTTl . the greatest amount of nonlinearity 
in the / — V characteristics is achieved with nearly perfect 
rotational orientation of graphene layers. This presents a 
significant challenge in fabrication of such devices based 
on the layer transfer technology. But epitaxial growth of 



graphene, BN |18| or other 2D materials provides a choice 
to overcome this problem. We note that due to the many- 
particle nature of the excitonic condensate, the BiSFET is 
expected to be insensitive to rotational misalignment [3|. This 
is similar to the robustness of superconductivity to defects. 
However, the single-particle tunneling nature also makes the 
SymFET robust to certain quantities to which the BiSFET 
is sensitive. As discussed earlier, the SymFET is robust to 
temperature. Another advantage is the robustness of single- 
particle tunneling - this is an intrinsic advantage for the 
SymFET. Although the tunneling current will vary with the 
tunneling insulator thickness and its dielectric constant, the 
regular single-particle tunneling behavior will survive in all 
temperature and thickness values. Thus, unlike the BiSFET's 
sensitivity to thickness variations, the SymFET behavior is 
robust to thickness and dielectric constant variations. Both 
thickness variations and inelastic scattering processes can be 
significantly suppressed by using 2D crystal 'insulators' such 
as BN or M0S2 between the graphene layers, preferably in a 
rotationally aligned structure. 

In summary, we have presented an analytical model to 
calculate the channel doping potential and current-voltage 
characteristics in a novel electronic device structure, the 
SymFET. The current in a SymFET flows by tunneling from 
one graphene layer to the other. The current is insensitive to 
temperature. The resonant current peak is controlled by chem- 
ical doping and applied gate bias. The on/off ratio increases 
with graphene coherence length and doping. The symmetric 
resonant peak is a good candidate for high-speed analog 
applications, of which frequency multiplication is an example. 
The resonant peak behavior can also be the framework for new 
digital architectures that consume much lower power than the 
current state of the art electronic switches. 
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